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1. Introduction 



The covariant quark model with infrared confinement developed in a series of papers (see 
Refs. [JXJ] — [J3J] ) is a successful tools for a unified description of the multiquark states: mesons, 
baryons, tetraquarks, etc. The covariant quark model is an effective quantum field approach to 
hadronic interactions based on the interaction Lagrangian between hadrons and their constituent 
quarks. Knowing a corresponding interpolating quark current allows calculating the matrix ele- 
ment of physical processes in a consistent way. A distinctive feature of this approach is that the 
multiquark states, such as baiyons (three quarks), tetraquarks (four quarks), etc., can be consid- 
ered and described as rigorously as the simplest quark-antiquark systems (mesons). The coupling 
constants between hadrons and their interpolating quark currents are determined from the connec- 
tion condition Z# = proposed in Refs. [|], ^ and used further in numerous subfields of particle 
physics (for a review, see Refs. [§, g Here Z H is a renormalization constant of the hadron 



wave function. The matrix elements of physical processes are determined by a set of associated 
quark diagrams, which are constructed according to 1 /N c — expansion. In the covariant quark model 
an infrared cutoff is effectively introduced in the space of Fock-Schwinger parameters, which are 
integrated out in the expressions for the matrix elements. Such a procedure allows one to eliminate 
all the threshold singularities associated with quark production and thereby ensures quark confine- 
ment. The model has no ultraviolet divergences due to vertex hadron-quark form factors, which 
describe a nonlocal structure of hadrons. The covariant quark model has a few free parameters: a 
mass of constituent quarks, an infrared cutoff parameter that characterizes confinement region, and 
parameters that describe an effective size of hadrons. 

We review here the last applications of the covariant quark model for studying the prop- 
erties of the B s — meson, the light baryons and tetraquarks. The form factors of the B(B S ) —> 
P(V)— transitions are evaluated in the full kinematical region of momentum transfer squared. As 
an application of the obtained results the widths of the B s — nonleptonic decays are calculated. The 
modes D~D+, D* D+ + DjD* + and D* D* + give the largest contribution to Ar for the B s - B s 
system. The mode J/\j/<p is suppressed by the color factor but it is interesting for the search of 
CP-violating New-Physics possible effects in the B s — B s mixing. 

The static properties of the proton and neutron, and the A-hyperon (magnetic moments and 
charge radii) and the behavior of the nucleon form factors at low momentum transfers are de- 
scribed. The conservation of gauge invariance of the electromagnetic transition matrix elements in 
the presence of a nonlocal coupling of the baryons to the three constituent quark fields is discussed. 

The consequences of treating the X(3872) meson as a tetraquark bound state are explored. The 
decay widths of the observed channels X — > J /y + 2n{'57i) and X — > D° +D + 7T° via the inter- 
mediate off-shell states X — > J /y + p(&>) and X — > D + D* are calculated. Its one-photon decay 
X — > y + J/yis also analyzed. The matrix element of the transition X — > y + J/yis calculated 
and its gauge invariance is proved. For reasonable values of the size parameter Ax of the X(3872) 
consistency with the available experimental data is found. The possible impact of the X(3872) in a 
s-channel dominance description of the J / y dissociation cross section is discussed. 
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2. Covariant quark model 

The coupling of a hadron H to its constituent quarks is described by the Lagrangian: 

£Lt = gH-H(x)-J H (x) (2.1) 
where the quark currents are defined as 



Jm{x) = jdx x j dx 2 F M (x,xi ,x 2 ) ■ q a fl (x\ ) T M q a fl (x 2 ) Meson 

Jb(x) = JdxiJdX2jdX3Fs(x,Xi,X2,X3) 



x T iq a f \{x x ) [q%(x 2 )Cr 2q %(x 3 ))-e a ^ Baryon 

Jt(x) = jdx\ .. . Jdx 4 F T (x,xi ,x 4 ) 

x (> (xi ) CTi q% (x 2 )) ■ (jc 3 ) T 2 Cctl M • £ fllfl2C £ fl3fl4C Tetraquark 



where T is a Dirac matrix or a string of Dirac matrices which projects onto the spin quantum 
number of the hadron H{x). The matrix C = y°y 2 is the usual charge conjugation matrix and the 
a, (i = 1,2,3) are color indices. The function Fh is related to the scalar part of the Bethe-Salpeter 
amplitude and characterizes the finite size of the hadron. To satisfy translational invariance the 
function Fh has to fulfil the identity Fh (x + a,x\ +a,... ,x n + a) = Fh(x,xi , . . . ,x n ) for any four- 
vector "a". In the following we use a specific form for the scalar vertex function 



F H (x,xi,...,x n ) = 8 yc-J^WjXjj ®H[Yj,(xi-Xj) 2 ) , (2.2) 
where <&h is the correlation function of the constituent quarks with masses m,-, (z = 1, . . . ,n) and 

n 

the mass ratios vv; = m,/ £ ntj. 

7=1 



The coupling constant gn in Eq. (2.1) is determined by the so-called compositeness condition 
originally proposed in Refs. [||, [7|] and extensively used in Refs. [[| ^, [K]]. The compositeness 
condition requires that the renormalization constant of the elementary hadron field H{x) is set to 
zero 

Z H = l-^n^(m|) = (2.3) 
where is the derivative of the hadron mass operator. To clarify the physical meaning of the 



compositeness condition in Eq. (|2.3|), we first want to remind the reader that the renormalization 



1/2 

constant Z ff can also interpreted as the matrix element between the physical and the corresponding 
bare state. The condition Z# = implies that the physical state does not contain the bare state 



and is appropriately described as a bound state. The interaction Lagrangian of Eq. (2.1) and the 
corresponding free parts of the Lagrangian describe both the constituents (quarks) and the physical 
particles (hadrons) which are viewed as the bound states of the quarks. As a result of the interaction, 
the physical particle is dressed, i.e. its mass and wave function have to be renormalized. The 
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condition Z# = also effectively excludes the constituent degrees of freedom from the space of 
physical states. It thereby guarantees that there is no double counting for the physical observable 
under consideration. The constituents exist only in virtual states. One of the corollaries of the 
compositeness condition is the absence of a direct interaction of the dressed charged particle with 
the electromagnetic field. Taking into account both the tree-level diagram and the diagrams with 
the self-energy insertions into the external legs (i.e. the tree-level diagram times Z# — 1) yields a 
common factor Zh which is equal to zero. 

We have used free fermion propagators for the quarks given by 

S t (k) = (2.4) 
nii— If 

with an effective constituent quark mass m,. 

For calculational convenience we will choose a simple Gaussian form for the vertex function 
^h{— k 2 ). The minus sign in the argument of this function is chosen to emphasize that we are 
working in Minkowski space. One has 

<t> H (-k 2 )=exp(k 2 /A 2 H ) (2.5) 

where the parameter A# characterizes the size of the hadron H. Since k 2 turns into — k\ in Eu- 



clidean space the form (2.5) has the appropriate fall-off behavior in the Euclidean region. We 
emphasize that any choice for <J># is appropriate as long as it falls off sufficiently fast in the ultra- 
violet region of Euclidean space to render the corresponding Feynman diagrams ultraviolet finite. 
As mentioned before we shall choose a Gaussian form for <J># for calculational convenience. 

We have included the confinement of quarks to our model in Ref. [[I]]. It was done, first, by 
introducing the scale integration in the space of a -parameters, and, second, by cutting this scale 
integration on the upper limit which corresponds to an infrared cutoff. In this manner one removes 
all possible thresholds presented in the initial quark diagram. The cutoff parameter is taken to 
be the same for all physical processes. We have adjusted other model parameters by fitting the 
calculated quantities of the basic physical processes to available experimental data. 

Let us give the basic features of the infrared confinement in our model. All physical matrix 
elements are described by the Feynman diagrams which are the convolution of the free quark 
propagators and vertex functions. Let n, I and m be the number of the propagators, loops and 
vertices, respectively. In Minkowski space the £-loop diagram will be represented as 

n m n 

U{p U ..., Pm ) = [d A k} ( n ®h+n (-4+n) 11 S h(k +Vfe), 

J i 1= l 1 3 =1 

4 + * = I(*&l + vg„) 2 (2.6) 
n 

where the vectors kj are linear combinations of the loop momenta k\. The v, are linear combinations 
of the external momenta p\ to be specified in the following. The strings of Dirac matrices appearing 
in the calculation need not concern us since they do not depend on the momenta. The external 

m 

momenta pi are all chosen to be ingoing such that one has £ p\ = 0. 

i=i 
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Using the Schwinger representation of the local quark propagator one has 

oo 

S(k) = (m+ jO / d P e ~ P {m2 ~ k2) (* 2 < ™ 2 ) ■ 
o 

For the vertex functions one takes the Gaussian form. One has 

& i+n ( -K 2 ) = exp \pi +n K 2 ] i = l,...,m, (2.7) 

where the parameters j3 1+ „ = Si = 1/A? are related to the size parameters. The integrand in Eq. ( [D5| ) 
has a Gaussian form with the exponential kak + 2kr + R where a is I x £ matrix depending on the 
parameter J3;, r is the ^-vector composed from the external momenta, and R is a quadratic form 
of the external momenta. Tensor loop integrals are calculated with the help of the differential 
representation 

0e 2kr = -^-e 2kr , 



We have written a FORM Qll[ ] program that achieves the necessary commutations of the differential 



operators in a very efficient way. After doing the loop integrations one obtains 

oo 

n = Jd n pF(p u ...,p n ), 



o 



where F stands for the whole structure of a given diagram. The set of Schwinger parameters j3, can 
be turned into a simplex by introducing an additional ^-integration via the identity 



f dt8{t-j^fr 

o ! ' =1 



leading to 



n = Jdtf- X Jd n a8(l-Y j a^F(ta u ...,ta n ). (2.8) 

o o 1=1 

There are altogether n numerical integrations: (n— 1) a-parameter integrations and the integration 
over the scale parameter t. The very large ^-region corresponds to the region where the singularities 
of the diagram with its local quark propagators start appearing. However, as described in [[I]], if one 
introduces an infrared cut-off on the upper limit of the t-integration, all singularities vanish because 
the integral is now convergent for any value of the set of kinematic variables. We cut off the upper 
integration at 1/A 2 and obtain 

l/A 2 1 

n c = J dtt n - l Jd n a8(\-Y i a^F{ta u ...,ta n ). 

i=l 

By introducing the infrared cut-off one has removed all potential thresholds in the quark loop 
diagram, i.e. the quarks are never on-shell and are thus effectively confined. We take the cut-off 
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parameter A to be the same in all physical processes. The numerical evaluations have been done 
by a numerical program written in the FORTRAN code. 

As a further illustration of the infrared confinement effect relevant to the applications in this 
paper we consider the case of a scalar one-loop two-point function. One has 

2 fd^E 

2{p > J K 2 [m 2 + {kE + y E) 2 ][m 2 + [kE _y E) 2 ] 

where we have collected all the nonlocal Gaussian vertex form factors in the numerator factor e~ sk % . 
Note that the momenta kg, p E are Euclidean momenta. Doing the loop integration one obtains 

i 



n 2 (p 2 



dt 



{s + tf 



da exp 



" fZloc + 



St 



Z] 



Zloc 



m 



a(l - a)p 2 



zi 



s+t 
1\ 2 



a ~2 



(2.9) 



The integral Yl2{p 2 ) can be seen to have a branch point at p 2 = Am 2 because zi oc is zero when 
a = 1/2. By introducing a cut-off on the f-integration one obtains 



n 2 ( P 2 



l/X 2 

dt 



o 



{s + tf 



da exp 



fZloc + ' 



St 

S + t' 



(2.10) 



The one-loop two-point function Yl c 2 (p ) Eq.(2.1C) can be seen to have no branch point at p =4m . 

The gauging of the nonlocal Lagrangian in Eq. (2. 1) proceeds in a way suggested in Refs. [ 12, 
13] and used before by us (see, for instance, Refs. [jl4, 15]). In order to guarantee local invariance 



of the nonlocal Lagrangian in Eq. (2T) one multiplies each quark field q(xi) with a gauge field 
exponential: 



qi{xi) 



,-ie q ,I(xj,x,P) 



qi{Xi) 



where 



I(xi,x,P) = J dz^iz). 



(2.11) 



(2.12) 



The path P connects the end-points of the path integral. One then expands the gauge exponential 
up to the requisite power of e q A^ needed in the perturbative series. We need to know only the 
derivatives of the path integral expressions when calculating the perturbative series. Therefore, we 
use the formalism suggested in [|12|, [13|] which is based on the path-independent definition of the 
derivative of I(x,y,P): 



J^J^TJ 1 ^^ = }ir n lI(x + dx,y,P')-I(x,y,P)} 



(2.13) 



where the path P' is obtained from P by shifting the end-point x by dx. The definition ( 2.13 ) leads 
to the key rule 



dxV 



I(x,y,P) =A fl (x) 



(2.14) 
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which in turn states that the derivative of the path integral I(x,y,P) does not depend on the path P 
originally used in the definition. 

As a result of this rule we are getting the part of the Lagrangian which describes the nonlocal 
interaction of the hadron, quark and electromagnetic fields to the first order in the electromagnetic 
charge. 



3. i^-meson 

We give below the necessary definitions of the leptonic decay constants, invariant form factors 
and helicity amplitudes. 

The leptonic decay constants of the pseudoscalar and vector mesons are defined by 

fpP^i P 2 = m 2 P , 

m v f v ei}, p 2 = ml, (3.1) 

where N c = 3 is the number of colors. 

Herein our primary subjects are the following matrix elements, which can be expressed via 
dimensionless form factors: 



N c gp 
N c gv 



d"k 
d 4 k 



o» P (-r)tr 



(2tt) 4 / 



O^Si{k + w i p)y 5 S 2 {k-W2p) 



( P to q3 )(P2)\42°' l 4l\ P l<h<lz]fa)) 

d 4 k 

r O^S l (k + p [ )y 5 S 3 {k)y 5 S 2 {k + p2 



Ncgpgp 



®p( -(k + Wl3Piy)<$>pi[-(k + W 2 3P2) 2 

--F + {q 2 )P»+F_(q 2 )q» 



(3.2) 



(P{g iq3] (P2)\q2(a flv q v )q l ^(pi)) 



Ncgpgp 



d 4 k 



®p[-{k + Wi3Pi) )®pi ( - (k + W23P2 



x tr 



a" V Si(&+/>i)y 5 Siifyy 5 S 2 (k + P2 ) 



(q 2 P^-q-Pq^)F T (q 2 ), (3.3) 



(V{p2,S2)[g iq3 ]\q 2 0^q l \P [ g m] (p l )) = 
Ncgpgv j j^^-. <&p ( - wi3j?i ) 2 ) ®v ( - (k + w 23 p 2 ) 2 ) 

r O" Si (k + Pl )fS 3 {k) ^S 2 {k + p 2 ) 



(3.4) 



fcy 

mi + ni2 



-g^P-qA (q 2 )+P^P v A + (q 2 )+q^P v A4q 2 )+ie^P a qpV(q 2 )), 
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(V(j>2,£2 
Ncgpgv 



<H13\ I 

d 4 k 



q2{a^q v {\ + ^))qx\P [q , q2] {px)) = 

- (^ + Wi3Pi) 2 ) 3>v( - {k + w 2 3P2) 2 



x tr 



(2x) 4 i 



(3.5) 



el -q^qVq 2 )P-qa,{q 2 ) + {P"P V -q"P V P-q/q 2 )a + {q 2 ) + ie^P a q^{q 2 )). 



Here, P = p\ + p2, q = p\ — pi, e| • p 2 = 0, = ra?. Since there are three sorts of quarks involved 
in these processes, we introduce the notation with two subscripts vv/j = m qj j(m qi +m q .) (i,j = 



1,2,3) so that Wjj + Wjj = 1. The form factors denned in Eq. ( p3| ) satisfy the physical requirement 
ao(0) = a+(0), which ensures that no kinematic singularity appears in the matrix element at q 2 = 0. 
For reference it is useful to relate the form factors we have denned to those used, e.g., in Ref. fll6|], 
which are denoted by a superscript c in the following formulae: 



f+, F 



mr — m% 



(/+ - fo ) > F T = fr: 



mi+W2 Af, A_ 
mi — m 2 



■A c 2 , A_ 



'2 > 



Tf, a + = T 2 c + 



2ni2(mi + m 2 ) 
q 2 



Ti. 



(A c 3 -A c ), V = V C , 



~ 2 2 J 3 
m l ~~ m 2 



(3.6) 



We note in addition that the form factors A\ (q 2 ) satisfy the constraints: Aq(0) = A c 3 (0) and 

2m 2 A c 3 (q 2 ) = [m x + m 2 )A\{q 2 ) -{m v - m 2 )A c 2 {q 2 ) . 

It is convenient to express all physical observables through the helicity form factors H m . The 
helicity form factors H m can be expressed in terms of the invariant form factors in the following 
way (see Refs. 0, [lSj, [l9|]): 



(a) Spin 5 = 0: 



Ho 



1 



vr 
o, 

2m! |p 2 | 



{(m 2 -m 2 )F + + q 2 F_} , 



(3.7) 



(b) Spin 5= 1: 



1 



mi |p 2 | 



mi +m 2 m 2 V^ 2 



{(m 2 -m 2 2 )(A + -A )+q 2 A_} 
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H± = \ {-(m?-mi)A ±2m 1 |p 2 |y}, (3.8) 

H Q = 1 — l - T ={-(m 2 l -ml)(m 2 l -ml-q 2 )A Q +4m 2 l \p 2 \ 2 A + } , 

where |p2| = A 1 1 2 (m^m^q 2 ) / \2m\) is the momentum of the outgoing particles in the rest frame 
of ingoing particle. 

The effective Hamiltonian describing the B s -nonleptonic decays is given by (see, Ref. [[!]]) 
v2 j = i 

Qi = (c ai b a2 ) V - A (s a2 c ai )v-A, Qi = {c a{ b ai ) V - A ,(s a2 c a2 )v-A, 

Q3 = (^i^aJv-A^CaJy-A, Q4 = (s ai b ai )v~A {c ai c ai )v-A , 

Qs = (s ai b ai )v-A(c a2 c a2 ) v+A , £>4 = {s ai b a2 ) v - A {c a2 c ai ) v +A , (3-9) 

where the subscript V — A refers to the usual left-chiral current = y^(l — y 5 ) and V +A to the 
usual right-chiral one 0+ = y^(l + y 5 ). The a, denote the color indices. 

We consider the nonleptonic decays of the S^-meson into D7 Z>+, D7 D* + , D* ~~ D+, D* ~ D* + 
and J/\j/<p. The calculation of the matrix elements is straightforward. It directly leads to the 
representation corresponding to naive factorization. 

The widths can be conveniently expressed in terms of the helicity form factors and leptonic 
decay constants. In the case of the color-allowed decays one has 

r(B s h. D-Dj) = \^ ] f { Cfm Ds M B?" (ml ) + 2Cfjg if* « 

16jt m\ 

r(a,->z>;-z>;+) = ||{^[A c ( ' v) ] 2 (cf m D . f D *f £ (flf^'Ko) 2 . O.io) 

Here, Ac = |V c foV^| and |q^| is the momentum of the second outgoing particle in the rest frame 
of B s — meson. The Wilson coefficients are combined as C| ff = C 2 + E, C\ + C4 + t, C3 and C| ff = 
C6 + % C5 . where a color factor = 1 /A^ will be suppressed in the numerical calculations according 
to 1 /N c — expansion. Also we do not take into account the annihilation channels which are available 
for the color-allowed decays. 

The width of the color-suppressed B s J/yQ decay is written as 

m^J/Yt) = ^^[&]^Cf + Cf) 2 {m J/w f I/r ) 2 £ (H?^(m 2 J/v )) 2 pM) 
io7t m Bs j=0± v / 

where the Wilson coefficients are combined as Cf ff = C\ + % C 2 + C3 + % C4 and C 5 eff = C5 + % C(,~ 



16k m\ 

l(B s ^D s D s ) = — -—Ac C 2 m D J Ds H t (m D )+2C 6 /n /> 5 m 

167T »ij V m £>; 

16 71 m B 

2 

£_ 

16jz 
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The first application of our relativistic quark model with infrared confinement to the descrip- 
tion of the physical observables was done in our paper [y. We have fitted the model parameters to 
the leptonic and radiative decay constants of both pseudoscalar and vector mesons. Then we have 
calculated transition form factors and the widths of the Dalitz decays and compared the results 
with available experimental data. Here we calculate the form factors describing the transitions of 
the heavy B{B S ) —mesons into the light ones, e.g. 71, K, p , K* , . These quantities are of great inter- 
est due to their applications to semileptonic, nonleptonic and rare decays of the B and B s — mesons. 
Basically, they are calculated within the light-cone sum rules (LCSR) in the region of large recoils 
(small transfer momentum squared). Our approach allows one to evaluate the form factors in the 
full kinematical regions including zero recoil. First, we update the model parameters by fitting them 
to the leptonic decay constants, see Table [TJ, and the widths of the radiative decays, see, Table ||. 
The results of the fit for the values of quark masses, the infrared cutoff and the size parameters are 
given in Eqs. ( 3.12| ), ( 3.13 ) and ( 3.14 ), respectively. 



m u m s m c nib A 



(3.12) 



0.235 0.424 2.16 5.09 0.181 GeV 



A n A K A D A Ds A B A Bs A Bc A p 



(3.13) 



0.87 1.04 1.47 1.57 1. 



1.95 2.42 0.61 GeV 



A a Aj/ V A K * Ad* A D * A b * A b * 
■ ■ (3-14) 

0.47 0.88 1.48 0.72 1.16 1.17 1.72 1.71 GeV 



In Figs. [T|f| we plot our calculated form factors in the entire kinematical region < q 2 < q 



2 



max - 



For comporison we also show the results obtained in the light-cone sum rules [26]. The figures 
highlight the wide range of phenomena accessible within our approach. 

As was suggested in Ref. p7[], one can check how well the form factors satisfy the low recoil 
relations among them. In Fig. || we plot the ratios 

o _W) R _T 2 (q 2 ) R _q 2 Hq 2 ) 

K l — ,w 7 , , K 2 — ■ i 9 ^ ; K 3 — —j ■ , 9 n • (3.1 J) 

which in the symmetry limit should be all of order 1 — (2a s / (3n) In (fi/rrib), i.e. near one. One can 
see that similar to the LCSR form factors, it works reasonably well for R\ and R 2 but not for Rt,. 

It is interesting to compare the behavior of the form factor calculated from the triangle loop- 
diagram with those from vector-dominance model (VDM). In the case of the B — n— transition, one 
has 

f VDM 9 I — ~ 2" 

mjp - q 1 
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Table 1: Leptonic decay constants fy (MeV) used in the least-squares fit for our model parameters. 



This work 



Other 



Ref. 



fx 


128.7 


130.4 ± 0.2 


f 

JK 


1 ^ 1 

1 JO. 1 


1 1 1 A Q 
1 JO. 1 ± U.O 


fo 


205.9 


206.7 ± 8.9 


fo s 


257.5 


257.5 ± 6.1 


fs 


191.1 


192.8 ± 9.9 


/b s 


234.9 


238.8 ±9.5 


/b c 


489.0 


489 ±5 


fp 


221.1 


221 ± 1 



This work Other Ref. 



10] 



fm 

u 

fj/y 
fx* 
fo* 
fo* 

fB* 



198.5 
228.2 
415.0 
213.7 
243.3 
272.0 
196.0 
229.0 



198 ± 2 [glj| 

227 ± 2 [gl]] 

415 ±7 [H) 

217 ±7 [HJ 

245 ± 20 [|5p 

272 ± 26 [||| 

196 ±44 [||| 

229 ±46 pi] 



Table 2: Electromagnetic decay widths (keV) used in the least-squares fit for our model parameters. 



Process 


This work 


Data 


n° ->■ 77 


5.06 x 10~ 3 


(7.7 ±0.4) x 10~ 3 


T] c ^77 


1.61 


1.8 ±0.8 


— > 7T ± 7 


76.0 


67 ±7 


ft) -»■ 7T°7 


672 


703 ± 25 


K *± _^ R ±y 


55.1 


50 ±5 


K*° ^K°j 


116 


116 ± 10 


D* ± ->D ± y 


1.22 


1.5 ±0.5 


J/v -»• n c y 


1.43 


1.58 ±0.37 



The curves are plotted in Fig. || One can see that they agree with quite good accuracy. That means 
the quark loop in some sense contains an information on the B*-pole. 

As an application of the obtained results we evaluate the widths of the B^-nonleptonic decays. 
The modes D~D+, D* D^ + D~D* + and D* D*+ give the largest contribution to Ar for the 
B s — B s system. The mode J/yf<j> is suppressed by the color factor but it is interesting for the search 
of CP-violating New-Physics possible effects in the B s — B s mixing. 
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B-7C: F + (q ) 



B-K: F + (q) 



10 15 20 25 

q'(GeV 2 ) 



B-Jt: F (q) 



B-K: F (q ) 



10 15 

q ! (GeV 2 ) 



Figure 1: Our results for the form factors appearing in Eqs. ( p.2[ ) & ( |3.3| ) - Left panel, B — K— transition; 
and right panel, B — K— transition. For comporison we plot the curves given by LCSR from Ref. 



For the CKM-matrix elements we use the values from [21] 



\V„d\ \Vus\ \Vub\ \v cd \ \v cs \ \v cb \ 



0.974 0.225 0.00389 0.230 0.975 0.0406 



(3.16) 



For the Wilson coefficients we take [28] 



Cx C 2 



C 3 



Q C5 



c 6 



(3.17) 



-0.257 1.009 -0.005 -0.078 0.000 0.001 



evaluated to next-to-next-to leading logarithmic accuracy in MS (NDR) renormalization scheme at 



the scale n =4.8 GeV [29]. 



We also need the values of the B s — 0— transition evaluated at q 2 = mj, . We give them in 



Table || and compare with results of Ref. [30]. 

Finally, we give our results for the branching ratios in Table |4[ One can see that there is good 
agreement with available experimental data. 
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B-p: A (q") 



B-p: V(q-) 




20 



B-p: A (q) 



B-p: T,(q-) 




20 



Figure 2: Our results for the form factors appearing in Eqs. (3.4) & (5.5) for B — p— transition. For 
comporison we plot the curves given by LCSR from Ref. [E6|. 



Table 3: The relevant B s — <j>— form factors at q 2 = mj/y calculated in our work. For comparison we give 
the results of Ref. [ph. 



This work Ref. 



A ^ m )/ W ) °- 37 0.42±0.06 
A i( m 7/v/) 048 0.38±0.06 



0.56 0.82±0.12 
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B-K*: KAsf) 



B-K*: V(q ) 




B-K*: A (ql 



B-K*: T,(q") 



20 



Figure 3: Our results for the form factors appearing in Eqs.(5.4) & (3.5) for B — A"*— transition. For 
comporison we plot the curves given by LCSR from Ref. [E^l. 



Table 4: Branching ratios (%) of the B s -nonleptonic decays calculated in our approach. 



Process 



This work 



B s -> DJD+ 
B S ^D-D* + +D*~D+ 
B S ->D* D* s + 
B s -tJ/yrQ 



1.65 
2.40 
3.18 
0.14 



Data [El 



1 04+ 029 



2.8 ±1.0 
3.1 ±1.4 
0.14 ±0.05 
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B s -*: A.(q) 



B s -0: V(q-) 




B -<|>: A (q) 



I 1 (GeV 2 ) 



Figure 4: Our results for the form factors appearing in Eqs. ( p.4[ ) & ( p.5| ) for B s — (^—transition. For 
comporison we plot the curves given by LCSR from Ref. [E6J. 



B-K*: the ratios R. (i= 1,2,3) 



LCSRB-K*: the ratios R. (i=l,2,3) 



13 14 



16 17 18 

q ! (GeV 2 ) 



13 14 



16 17 18 19 

q" (GeV 2 ) 



Figure 5: Our results for the ratios of the form factors appearing in Eq. (3.15) for B — K*— transition 
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B-7t form factor from loop and B -monopole 



loop 

monopole 




I , 1 1 1 , 1 , 1 1 1 , 1 

5 10 15 20 25 30 

q 2 (GeV 2 ) 

Figure 6: The comparison of the results for the B — n— form factor obtained on the one hand from the 
quark-loop diagram and on the another hand from the VDM-monopole. 
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4. Light baryons 

Let us begin our discussion with the proton. The coupling of a proton to its constituent quarks 
is described by the Lagrangian 

^ t (x)=g N p(x)-J p (x)+g N J p (x)-p(x), (4.1) 



(4.2) 



where we make use of the same interpolating three-quark current J p (J p ) as in Ref. [14] 

J p (x) = JdxiJdx 2 Jdx 3 F N (x;x u x 2 ,x 3 )J^ i > (x l ,x 2 ,x 3 ), 

4*\x h x 2 ,x 3 ) = r A ^d^(x l )-[e a ^u^(x 2 )cr A u^(x 3 )], 

Jp{x) = j dx x J dx 2 J dx 3 F N {x\x u x 2 ,x 3 )j\ p ^ (x l ,x 2 ,x 3 ), 

J%\x u x2,x 3 ) = [£ fl ' fl2a3 «- fl3 (^)r A c« fl2 (^)]-^ 1 ^i)7 5 r A . 

The matrix C = y°y 2 is the usual charge conjugation matrix and the a\ (i = 1 , 2, 3) are color indices. 
There are two possible kinds of nonderivative three-quark currents: P 4 <g> Fa = Y* ® 7a (vector cur- 
rent) and r A (g>r A = j o a ^ ® o a n (tensor current) with o a ^ = ^(f^yP — yPy"). The interpolating 
current of the neutron and the corresponding Lagrangian are obtained from the proton case via 
p — > n and As will become apparent later on, one has to consider a general linear superpo- 

sition of the vector and tensor currents according to 

J N =x4 + (l-x)4, N = p,n (4.3) 

The electromagnetic vertex function A p (p, p') of the proton consists of four pieces represented 
by the four two-loop quark diagrams in Fig. 

Let us briefly describe a check on the gauge invariance of our calculation. Without gauge in- 
variance there are three independent Lorentz structures in the electromagnetic proton vertex which 
can be chosen to be 

^ P (p,p') = fF}\q 2 )-^-Fi(q 2 )+q^F 1 P G (q 2 ), (4.4) 

where o M = ^{y^y v — y v y^)qv The form factor F^ G (q 2 ) characterizes the non-gauge invariant 
piece and must therefore vanish for any q in a calculation which respects gauge invariance. For 
the four contributions of Fig. 2a-2d we found that 

^ Gd (<7 2 )=0, F£ Gu (q 2 )^0, ^ G(i ,)(? 2 ) = -^ G(fl) (? 2 ) V- (4.5) 

It means that the non-gauge invariant contributions of the two vertex diagrams are zero while they 
vanish for the sum of the two bubble diagrams. 

The electromagnetic vertex function of the neutron is obtained from that of the proton by 
replacing m u <-> mj, e„ f-> e c i and m p — > m n . F^(q 2 ) and F 2 N {q 2 ) are the Dirac and Pauli nucleon 
form factors which are normalized to the electric charge and anomalous magnetic moment 
(A:jv is given in units of the nuclear magneton e/2m p ), respectively, i.e. one has F^(0) = and 
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Figure 7: Electromagnetic vertex function of the proton: (a) vertex diagram with the e.m. current attached 
to d-quark; (b) vertex diagram with the e.m. current attached to u-quark; (c) bubble diagram with the e.m. 
current attached to the initial state vertex; (d) the bubble diagram with e.m. current attached to the final state 
vertex. 

i*2 (0) = kf]. In particular, one can analytically check by using the integration-by-part identity that 
the Dirac form factor of the neutron is equal to zero at q 2 = 0. 

The nucleon magnetic moments jx N = F^(0) + F 2 A '(0) are known experimentally with high 



accuracy [21] 



^ = 2.79 jU„ expt = - 1 -91 . (4.6) 
We will use these values to fit the value of the nucleon size parameter. We obtain 

vector current A N = 0.36GeV \i p = 2.79 \i n = -1.70, (4.7) 

tensor current =>• A N = 0.61GeV \i p = 2.79 \l n = -1.69. (4.8) 

It is convenient to introduce the Sachs electromagnetic form factors of nucleons 

G N E (q 2 ) = W) + ^/"(i) , G N M (q 2 ) = W) + W) • (4-9) 
The slopes of these form factors are related to the well-known electromagnetic radii of nucleons: 



2 \N 



-dG E N (q 2 " 



(4) 



dq 2 



2 \N 



(Hi) 



6 dG N M {q 2 ) 



^ — G^(0) dq 2 



(4.10) 

q 2 =0 



We would like to emphasize that reproducing data on the neutron charge radius (r|)" is a 



nontrivial task (see e.g. discussion in Ref.[31]). As well-known the naive nonrelativistic quark 
model based on SU(6) spin-flavor symmetry implies = 0. The dynamical breaking of the 
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SU(6) symmetry based on the inclusion of the quark spin-spin interaction generates a nonvanishing 
value of {r\) n . From this point of view the dominant contribution to the (r§)" comes from the Pauli 
term: 

<^4 Fi ' (o> ' 

The experimental data on the nucleon Sachs form factors in the space-like region Q 2 = —q 2 > 
can be approximately described by the dipole approximation 



E l+fl,, ft, « 2 ft, (l-,V0.71GeV 2 ) 2 



According to present data the dipole approximation works well up to 1 GeV 2 (with an accuracy 
of up to 25%). For higher values of Q 2 the deviation of the nucleon form factors from the dipole 
approximation becomes more pronounced. In particular, the best description of magnetic moments, 
electromagnetic radii and form factors is achieved when we consider a superposition of the V- and 



T-cunents of nucleons according to Eq. ( \[.3\ ) with x = 0.8. For the size parameter of the nucleon 
we take A N = 0.5 GeV. 

In Table || we present the results for the magnetic moments and electromagnetic radii for this 
set of model parameters. In Fig. || we present our results for the q 2 dependence of electromagnetic 
form factors in the region Q 2 £ [0, l]GeV 2 . Fig. ^ also shows plots of the dipole approximation 
to the form factors. The agreement of our results with the dipole approximation is satisfactory. 
Inclusion of chiral corrections as, for example, developed and discussed in [[32|] may lead to a 
further improvement in the low Q 2 description. 



Table 5: Elecnomagnetic properties of nucleons. 



Quantity 


Our results 


Data [21] 


\l p (in n.m.) 


2.96 


2.793 


\l n (in n.m.) 


-1.83 


-1.913 


rl (fm) 


0.805 


0.8768 ± 0.0069 


{r 2 ) n (fm 2 ) 


-0.121 


-0.1161 ±0.0022 


>M (fm) 


0.688 


0.777 ± 0.013 ± 0.010 


^(fm) 


0.685 


a o to +0.009 
U.SO2_ 008 
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G (mixing) 




G M / u (mixing) 



(4m / q") (G / a ), q= - Q", (mixing) 




G (mixing) 




Figure 8: Sachs nucleon form factors in comparions with the dipole representation in the space-like region 
Q< 1 GeV 2 . 



20 



Multiquark states in the covariant quark confinement model 



Mikhail A. Ivanov 



5. The X(3872) -meson as a tetraquark 

A narrow charmonium-like state X (3872) was observed in 2003 in the exclusive decay process 
B ± ->■ K+n+n'J /y Q. TheX(3872) decays into n + %-J/y and has a mass of m x = 3872.0 ± 
0.6(stat) ±0.5(syst) very close to the M D o +M D ,a = 3871.81 ±0.25 mass threshold [J2TJ]. Its width 
was found to be less than 2.3 MeV at 90% confidence level. The state was confirmed in B-decays 



by the BaBar experiment [34] and in pp production by the Tevatron experiments [35]. 

From the observation of the decay X(3872) — > J / yy reported by [[36|], it was shown that the 
only quantum numbers compatible with the data are J PC = 1 ++ or 2 v . However, the observation 
of the decays into D°D 7t° by the Belle and BaBar collaborations allows one to exclude the 
choice 2~ + because the near-threshold decay X — > D°D°^: is expected to be strongly suppressed 
for J = 2. 

The Belle collaboration has reported evidence for the decay mode X — > 7t + 7t~7t Q J /y with a 



strong three-pion peak between 750 MeV and the kinematic limit of 775 MeV []36|], suggesting that 
the process is dominated by the sub-threshold decay X — > coJ/y. It was found that the branching 
ratio of this mode is almost the same as that of the mode X — > 7i + n~J/y: 

^V J/ r VK+ T K ? = 1.0±0.4 ( stat)±0.3 ( syst). (5.1) 

These observations imply strong isospin violation because the three-pion decay proceeds via an 
intermediate co-meson with isospin whereas the two-pion decay proceeds via the intermediate 
p-meson with isospin 1. Also the two-pion decay via the intermediate p-meson is very difficult to 
explain by using an interpretation of the X(3872) as a simple cc charmonium state with isospin 0. 

There are several different interpretations of the X(3872) in the literature: a molecule bound 
state (D°D*°) with small binding energy, a tetraquark state composed of a diquark and antidiquark, 
threshold cusps, hybrids and glueballs. A description of the current theoretical and experimental 



situation for the new charmonium states may be found in the reviews [38]. 



We provided in Ref. [||] an independent analysis of the properties of the X (3872) meson which 




we inteipret as a tetraquark state as in [[39]]. The authors of [ p9| ] suggested to consider the X(3872) 
meson as a J PC = 1 ++ tetraquark state with a symmetric spin distribution: [cg].s = o [c#]s = i + [cq]s=\ [cq]s=o, 
(q = u,d). The nonlocal version of the four-quark interpolating current reads 

J xS x ) = Jdx x . 

x -^=S abc S dec | [q a (x 4 )Cfc b (x l )}[q d (x 3 )Y >1 Cc e (x2)} + (f O J* 1 ) }, (5.2) 

where w\ = W2 = m c /2(m q + m c ) and W3 = W4 = m q /2(m q +m c ). The matrix C = y°y 2 is the charge 
conjugation matrix. The effective interaction Lagrangian describing the coupling of the meson X q 
to its constituent quarks is written in the form 

^im = gx Xqn(x) ■ Jx q (x) , (q = u,d). (5.3) 

The state X u breaks isospin symmetry maximally so the authors of [|9|] take the physical states to 
be a linear superposition of the X u and Xj states according to 

X/ = Xi ow = X„ cos G + X c i sin G , 
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Xh = ^high = — X„ sin Q +Xd cos 9- 



(5.4) 



The mixing angle 6 can be determined from fitting the ratio of branching ratios Eq. ( p.l| ). 



The coupling constant g x in Eq. ( p.3| ) will be determined from the compositeness condition: 

Z x = l-W x (m 2 x )=0, 

where Yl x (p 2 ) is the scalar part of the vector-meson mass operator. The corresponding three-loop 
diagram describing the X-meson mass operator is shown in Fig. |[ 




u 




Figure 9: Diagram describing the X„-meson mass operator. 



Next we evaluate the matrix elements of the transitions X — >• J/y + p((o) and X — > D + D*. 
The relevant Feynman diagrams are shown in Fig. [To[ 




Figure 10: Feynman diagrams describing the decays X — > J/\j/+p((o) andX — > D + D*. 



Since the X(3872) is very close to the respective thresholds in both cases, the intermediate p, 
CO and D* mesons have to be treated as off-shell particles. Using the calculated matrix elements for 
the decay X — > J + p(oo) one can evaluate the decay widths X — > J/Y + 27i(3n). We employ 
the narrow width approximation for this purpose. 



There are two new free parameters: the mixing angle 6 in Eq. ( |5.4| ) and the size parameter A x . 
We have varied the parameter Ax in a large interval and found that the ratio 



r(X u ^J/y + 3 7i) 
T{X u ^J/y + 2n) 



0.25 



is very stable under variations of A x . Hence, by using this result and the central value of the 
experimental data given in Eq. (5.14), one finds 6 « ±18.4° for X/ ("+") and X/, ("-"), respectively. 
This is in agreement with the results obtained in both [39]: 6 w ±20° and [40]: 6 w ±23.5°. The 
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decay width is quite sensitive to the change of the size parameter Ax- A natural choice is to take a 
value close to Aj/^ and A^ c which are both around 3 GeV. We have varied the size parameter Ax 
from 2.4 up to 4 GeV and found that the decay width F(X — > J / y+ n 7i) decreases from 0.25 MeV 
monotonously. This result is in accordance with the experimental bound T(X(3872)) < 2.3 MeV 



and the result obtained in [J39[] : 1.6 MeV. 

In a similar way we calculate the width of the decay X — > D°D 7r° which was observed by the 
Belle Coll. and reported in [f37[]. As in the previous case we have varied Ax from 2.5 up to 4 GeV 
and found that the decay width T(X/ — > D°D 7t°) decreases from 1.1 MeV monotonously. We plot 



the dependence of the calculated decay widths on the size parameter Ax in Fig. 1 1 



r(X -> I/y + nil), MeV 




A v (GeV) 



Figure 11: The dependence of the decay widths T(Xi -> D°D Q 7t°) and T(X -> J /y + nn) on the size 
parameter Ax. 



Using the results of [21], one calculates the experimental rate ratio 

T(X -> D°D 7i ) 

—± ; —A- = 10.5±4.7 

r(X4J/p+r) 

The theoretical value for this rate ratio depends only weakly on the size parameter Ax : 



(5.5) 



T(X -> D°D°n Q ) 



T(X^J/xi/k + 7I- 



4.5 ±0.2. 



theor 



(5.6) 



The theoretical error reflects the Ax dependence of the ratio. The ratio lies within the experimental 



uncertainties given by Eq. (5.5). 



The matrix element of the decay X(3872) — ^ J/y + j can be calculated from the Feynman 



diagrams shown in Fig. [12. 

The invariant matrix element for the decay is given by 



M(X q (p) ^J/ W ( qi ) + y( q2 )) =i(2n) A d^\p-q l -q 1 )e^e p y e v J/xi ,T llpv {q u q 2 ) (5.7) 
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J/V 





Figure 12 : Feynman diagrams describing the decay X — > J / y/ + y. 

We have analytically checked on the gauge invariance of the unintegrated transition matrix 
element by contraction with the photon momentum q 2 which yields q^T tipv {qi,q 2 ) = using the 
identities 

S(k 2 ) 4iS(k 2 + q 2 ) =S(k 2 + q 2 )-S{k 2 ), 
l 

J dT&(-za - (1 - x)b) (a-b)= *(-&) - 
o 

If one takes the on-mass shell conditions 

e|p M =0, ej /v q lv =0, e$q 2p =0 (5.8) 
into account one can write down five seemingly independent Lorentz structures 

Tfip V {q\,q 2 ) = £q 2 nvpq 2 \ W\ +£ qm vpqiiiW 2 + £ qm p_ p q 2v W3 + £ qiq2 p. v q\p W4 + S q] p_ vp qiq 2 W5 . 
Further, using the gauge invariance condition 

q 2 T npv = qiq2£ qiq2 nv(W4 + w 5 ) = o 

one has W4 = — W5 which reduces the set of independent covariants to four: 

Tjipv{q\-,q 2 ) = £q 2 nvpqiW\ + £ qiq2 vpq\p. W 2 + £ qiq2 n P q 2v W3 + {e qiq2 p. v q\p — qiq 2 £ qi p.vp . 

The gauge invariance condition W4 = — W5 provides for a numerical check on the gauge invariance 
of our calculation as described further on. 
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However, there are two nontrivial relations among the four covariants which can be derived by 
noting [||] that the tensor 



7 )l[VlV2V 3 V4V 5 ] =§/iV 1 £v 2 V3V4V 5 + cvcl. ( Vi V 2 V 3 V 4 V 5 ) 



(5.9) 



vanishes in four dimensions since it is totally antisymmetric in the five indices (Vi, V 2 , V3, V4, V5). 
Upon contraction with q^q^q^ 2 an d cfyq^q^ 2 one finds 

qi£q 2 nvp + £q\q2vpq\\i + (zqiqinvqip — qiqiZqi 

q\qi£qi\ivp ~ £q\q2Vpq\jl ~ £q\qip.pq2v = • 

It reduces the set of independent covariants to two. This is the appropriate number of independent 
covariants since the photon transition is described by two independent amplitudes as e.g. by the E 1 
and M2 transition amplitudes. 

The quantities W; are represented by the four-fold integrals 



(5.10) 



o o 



where we have suppressed the additional dependence of the integrand F* on the set of variables 



p 2 ,q\,q\;m q ,m c ,sx,sj/y with sx = 1/A| and Sj/y = l/A 2 -,^. The integrals in Eq. ( 5.10 ) have 



branch points at p 1 = 4(m q +m c ) 2 (diagram in Fig. [l2|-a) and at p 2 = Am 2 (diagrams in Figs. 12- 
b,c,d). At these points the integrals become divergent in the convential sense when t — > oo. Under 
numerical check on gauge invariance of the amplitude Tup V (qi,qz), we assume that the X-meson 
momentum squared is below the nearest unitarity threshold, i.e. p 2 < 4m 2 . The gauge invariance 
condition is independent of the overall couplings gx and gj/ w and thus the numerical check can be 
done irrelevant of their values. 

In the next step we introduce an infrared cutoff 1/A 2 on the upper limit of the t-integration 



in Eq. (5.10). In this manner one removes all possible singularities and thereby guarantees quark 



confinement. However, the contributions coming from the bubble diagrams in Figs. 12-b,c,d blow 
up at p 2 = m\ compare with the contribution from the diagram in Fig. 12 a. The bubble diagrams 
are needed only to guarantee the gauge invariance of the matrix element. For physical applications 
one should take into account only the gauge invariant part of the diagram in Fig. 12 -a. 



It is convenient to present the decay width via helicity or multipole amplitudes. One has 



r(X^J/ ¥+Y ) 



1 |Q2 1 



\2% 



m\ 



HA +|//r| 2 



1 I Q2 1 



\2% 



where the helicity amplitudes Hi and Hj are expressed in terms of the Lorentz amplitudes as 

m x 



|A £ i| 2 + |A M 2| 



(5.11) 



H L = imxmj/ Y \q 2 \ 



Wi + - 



H T 



-im 



J/V 



02 1 



Wi + 



m Jlw 
m x 



<l2 W3 — W4 



m 



q 2 \W 2 -[l + 



J/V 



mxml 



m 



J/V 



W 4 



ml 



mi- 



ni 



J/V 



2mx 



(5.12) 
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r<X -> JA|f + 2n), McV 
r(X->J/\|l + Y), McV 



A.. (GeV) 



Figure 13: The dependence of the decay widths F(X[ 
parameter Ax. 



J/W+j) and r(X/ J /y + 2n) on the size 



Proceeding in a such way one will get the dependence of the decay widths r(X/ — > J/y+ /) an d 
r(X/ -> J/Y + 2n) plotted in Fig. ||. 

Note that the radiative decay width for Xh = —X u sin d +Xj cos is almost an order of mag- 
nitude less than for Xi = X u cosd +Xdsind. If one takes Ax € (3,4) GeV with the middle point 
Ax = 3.5 GeV then the ratio of the widths is equal to 



r(X,^J/y+Y) 



0.15 ±0.03 



theor 



(5.13) 



F(Xi^ J/y + 2n) 

which fits very well the experimental data from the BELLE collaboration written down in their 
Eq. ( gig ). 



TjX^J/y + y) 
T(X -^J/y+ln) 



0.14 ±0.05 BELLE tffijfr 
0.22 ±0.06 BABAR[i2| 



(5.14) 



The last topic which we would like to discuss is the impact of the intermediate X-resonance 
on the value of the //i//-dissociation cross section, see [ff3]|-[H]. The relevant s-channel diagram 
is shown in Fig. [14]. 

We take Tx = 1 MeV in the Breit-Wigner propagator and set A^ = 3.5 GeV when calculating 



the matrix elements. We plot the behavior of the relevant cross sections in Fig. [15J. One can see that 
in the case of charged D-mesons (left panel in Fig. |l5|) the maximum value of the cross section is 
about 0.32 mb at E = 3.88 GeV. This result should be compared with the result of the cross section 
a(J/^+n ->Z) + D*)«0.9mbat£ = 4.0 GeV, see, [||] and the result of the cross section 
a (J/y + p -)• D + D*) w 2.9 mb at E = 3.9 GeV, see, [||]. Thus the X-resonance gives a sizable 
contribution to the // ^-dissociation cross section. It would be interesting to do a complete analysis 
of the J/y dissociation cross section in view of our new results on the s-channel contribution of 
the X(3872) tetraquark state. 
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p(w) 




Figure 14: Diagram describing the X-resonance contribution to the //(//-dissociation process. 



Charged D-mesons 



Neutral D-mesons 





Figure 15: The cross sections of the processes J/y/ + v° — >• X — >• D + D* . Charged D-mesons- left panel, 
neutral D-mesons-right panel. 
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6. Summary 

• We have presented a refined covariant quark model which includes infrared confinement of 
quarks. 

• We have calculated the transition form factors of the heavy B s — meson to light pseudoscalar 
and vector mesons, which are needed as ingredients for the calculation of the semileptonic, 
nonleptonic, and rare decays. Our form factor results hold in the full kinematical range of 
momentum transfer. 

• We have made use of the calculated form factors to calculate the nonleptonic decays B s — > 
D S D S ,... and B s — > J/ycj), which have been widely discussed recently in the context of B s — 
B 4 -mixing and CP violation. 

• We have applied our approach to baryon physics by using the same values of the constituent 
quark masses and infrared cutoff as in meson sector. 

• We have calculated the nucleon magnetic moments and charge radii and also electromagnetic 
form factors at low energies. 

• The properties of the X(3872) as a tetraquark have been studied in the framework of a co- 
variant quark model with infrared confinement. 

• The matrix elements of the off-shell transitions X — > J /y + p{co) and X — > D + D* were 
calculated. 

• The obtained results were then used to evaluate the widths of the experimentally observed 
decays X->J/y + 2k(3k) andX ^ D Q + D° + n°. 

• The possible impact of the X (3872) on the //^-dissociation process was disscussed. 

• We have calculated the matrix element of the transition X — > y + J/y and have shown its 
gauge invariance. We have evaluated the X — > y+J/y decay width and the polarization of 
the // Y in the decay. 

• The comparison with available experimental data allows one to conclude that the X(3872) 
can be a tetraquark state. 
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